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ABSTRACT 

We present the first three-dimensional simulations of fingering convection per- 
formed in a parameter regime close to the one relevant for astrophysics, and 
reveal the existence of simple asymptotic scaling laws for turbulent heat and 
compositional transport. These laws can straightforwardly be extrapolated to 
the true astrophysical regime. Our investigation also indicates that thermocom- 
positional "staircases," a key consequence of fingering convection in the ocean, 
cannot form spontaneously in the fingering regime in stellar interiors. Our pro- 
posed empirically-determined transport laws thus provide simple prescriptions 
for mixing by fingering convection in a variety of astrophysical situations, and 
should, from here on, be used preferentially over older and less accurate pa- 
rameterizations. They also establish that fingering convection does not provide 
sufficient extra mixing to explain observed chemical abundances in RGB stars. 

Subject headings: hydrodynamics — instabilities — turbulence 



1. Introduction 

Recent years have seen rapidly growing interest in astrophysical applications of fingering 
convection (often referred to as "thermohaline convection"), a process by which the unequal 
diffusion rates of thermal and compositional fields can drive a double-diffusive instability 
and lead to significant turbulent transport across stably-stratified regions. In the Earth's 



-2- 



oceans th is process ha s been well-studied as "salt fingering" (see ISternl Il960l . or a recent 



review m lKunzdl2003l ). In stellar contexts, by contrast, its exact form and contribution to 
vertical mixing is much less clear. Nevertheless, fingering convection has been invoked in 
various scenarios. It is thought to explai n the lack of observed metal licity signatures from 
the infall of planets onto their host star ( jVauclairl 12004 ; iGaraudl 120101 ) . It may also explain 
the presence of extra mixing in the radiative zones of RGB stars needed to fit abundance 
observations (ICharbonnel fc Zahnl 120071 ; IStancliffd 12010 ; IDenissenkovl |2010| ) . 



While field, laboratory, and numerical measurements abound for salt fingering, data 
is much more scarce for the astrophysical regime. Transport by fingering convection in 



astrophysics has until now bee n mod e 



ed th rou gh mixing-length theory, commonly used 



parameterizations being those of lUlrichl (119721 ) and lKippenhahn et al.l (ll980l ). In these works 



turbulent compositiona l trans port depends sensitively on the finger aspect ratio. As recently 
shown by IDenissenkovl (120101 ). however, there is a fundamental contradiction between the 
large aspect r atio required to explain a bundance observations of RGB stars using these 
prescriptions (ICharbonnel fc Zahnl 120071 ) . and the near- unit aspect ratio of the "fingers" 
actually measured in his own two-dimensional simulations. 



Three possib ilities emerge to resolve this problem. First, the iKippenhahn et al.l (ll980l ) 
and lUlrichl ( 119721 ) parameterizations may not adequately model transport by fingering con- 
vection. Secondly, two-dimensional simulations may not correctly capture the true dynamics 
of this inherently three-dimensional system. Finally, thermocompositional layers may form in 
the astrophysical objects studied, in which case transport could be strongly enhanced above 
the level of homogeneous fingering convection. Indeed, "staircases" of convectively mixed 
layers separated by thin fingering interfaces are ubiquitous in the ocean. In their presence, 
vertical mixing increases by as much as an o r der of magnit ude above the already-enhanced 
mixing due to fingering (jSchmitt et al.l 120051 ; IVeronisI 120071 ). Whether such staircases form 
in astrophysics has never been established. 

In this Letter we consider all three possibilities and present the first three-dimensional 
simulations to address directly the question of thermal and compositional transport by fin- 
gering convection in astrophysics. Using a high-performance algorithm specifically designed 
for homogeneous fingering convection, we are able to run simulations at moderately low 
values of the Prandtl number (Pr = v/kt) and diffusivity ratio (r = k^/kt), where v is 
viscosity, and kt and are the thermal and compositional diffusivities respectively. We 
find, in Section I3.1[ asymptotic scaling laws for transport which are applicable to the stellar 
parameter regime (Pr <C l,r <C 1). These provide a parameter- free, empirically- motivated 
prescription for mixing by homoge neous fingering co nvection. The compositional turbulent 
diffusivity derived is, as expected (jStern et al.ll200ll ) a few times larger than that obtained 
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by iDenissenkovl (120101 ) for the same parameters, but remains too small to account for RGB 
abundance observations. 

Since the presence of thermocompositional staircases is known to increase transport 
in the ocean, we next address the question of whether such staircases are likely to form 
spontaneously in the astrophys ical context, in Section 13.21 Apply ing a recently-validated 
theory for staircase formation (jRadkol 120031 : IStellmach et al.ll2010l ). we find that they are 
in fact not expected to appear in this case. We conclude in Section @] that our small-scale 
flux laws can reliably be used to estimate heat and compositional transport for a variety of 
astrophysical fingering systems and concur with Denissenkov's view that fingering convection 
is not sufficient to explain the required extra mixing in RGB stars. 



2. Model description and numerical algorithm 



The dynamics of double-diffusive systems depend on the fluid properties (u, k t , k^), as 
well as its local stratification, measured by V — V a d (where V = d\iaT/d\np and V a d = 
(d\iaT/d\np) ad is the adiabatic gradient) and (where = dln/x/dlnp). Here T is 
temperature, p is pressure, and \i is the mean molecular weight. The fingering instability 
occurs in thermally stably-stratified systems (V — V a d < 0) destabilized by an adverse 
compositional gradient (V M < 0). In most regi mes of intere st the main governing parameter 
is the density ratio, defined in astrophysics as (jUlrichlll972l ): 

V- V ad 



tin 



(1) 



This l inear instability is w ell understood, thanks to early work i n the o c eanic conte xt (IStern 



960; iBaines &: Gilll Il969[ ). later extended in astrophysics by lUlrichl (119721 ) and ISchmitt 



(119831 ) . A system is Ledoux- unstable when R$ < 1, and fingering-unstable when Rq e 
[1, 1/t]. Unfortunately, as discussed in Section [TJ the nonlinear saturation of this instability, 
or in other words its turbulent properties, remain poorly understood. 

We approach the problem from a numerical point of view, and model a fingering-unstable 
region using a local Cartesian frame (x,y,z) with gravity g = —ge z . We run our numeri- 
cal experiments using a high-performance spectral code, which was recently used to model 
three-dimensional fingering convecti on at oceanic parame ters (ITraxler et al.ll2010l ). as well 
as thermohaline staircase formation (IStellmach et al.ll2010l ). As a consequence of its original 
purpose, our code uses the Boussinesq approximation, which relates the density perturba- 
tions p to the temperature and compositional perturbations T and fx via 



— = —of + (3u, 
Po 



(2) 
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where po is a reference density, a = —p^dp/dT, and (3 = p$ dp/dp,. Fingering convection is 
driven by constant large-scale temperature and compositional gradients T 0z and p 0z . These 
approximations, which greatly simplify the governing equations, are nevertheless justified 
since fingers are typically much smaller than a pressure or densi ty scale heigh t, and velocities 
are small compared with the sound speed. As first shown by lUlrichl (119721 ). results in the 
Boussinesq case can straightforwardly be extended to the astrophysical case by replacing the 
Boussinesq density ratio R = aT 0z / (3p 0z with the true density ratio -Rq defined in equation 

We now describe the governing equations in more detail. We express the velocity, tem- 
perature, and composition fields as the linear background stratification plus perturbations, 



T(x,y,z,t) 
p(x,y,z,t) 
u(x,y,z,t) 



T (z) +f(x,y,z,t), 
p (z) + p(x,y,z,t), 
u(x,y,z,t), 



(3) 
(4) 
(5) 



where T (z) = T 0z z and po(z) = PozZ- We nondimensionalize the system noti ng that an 
appropriate length scale is provided by the anticipated finger width ( ]Sternlll960l ). [/] = d = 
{k t v / gaTQz) 1 / 4 . We then scale time using d 2 /n T , temperature with T 0z d, and composition 
with (a/ (3)T 0z d. The nondimensional governing equations are: 



1 fdu 



Pr \dt 
df 



+ u ■ Vu 
Vu 

+ w + u • VT 



dt 

dp 1 _ _ 



-Vp + (f - p)e z + V 2 u, 
0. 

v 2 f, 



rV 2 /i, 



(6) 
(7) 
(8) 

(9) 



Note that the Rayleigh number Ra then depends only on the computational domain height 
L z : 

®gTo z L 4 z ( L, 



Ra 



d 



We solve (EH9J) in a triply-periodic box of size (L x , L y ,L z ), so that 

f(x, y, z, t) = f(x + L x , y, z, t) = f(x, y + L y , z, t) = f(x, y } z + L z , t), 
and similarly for p, p and u. 



(10) 



1 Because the length scale of the convective motions is set by the diffusive scales in the fingering problem, 
it does not suffer from th e known pathology of t riply-periodic thermal convection, i.e. the homogeneous 
Rayleigh-Benard problem ( Calzavarini et al. 20061 ). 
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3. Results 

3.1. Small-scale fluxes 

Six sets of simulations were conducted at successively lower values of Pr and r, for a 
range of density ratio values l isted in Table [TJ spa nning in each case the instability range 



Ro G [1,1/r]. As shown by iTraxler et al.l (120101 ) . the statistical properties of fingering 



convection are independent of the computational domain size (i.e. of Ra and aspect ratio) 
provided it is wide enough to contain a sufficient number of fingers. The domain must also be 
tall enough to ensure that individual fingers do not span the entire height of the box. In most 
cases, the simulation domains were 67c? x 67 d x 107.2c? with a resolution of 96 3 grid points. 
At the selected parameter values, this correspo nds to approximately 5 x 5 x 10 wavelengths 



of the fastest-growing mode (see ISchmittl 119791 ). thus allowing space for many fingers to fit 



in the domain. As Rq — > 1/r, where fingers tend to become more vertically elongated, wider 
and taller simulation domains (83.75c? x 83.75c? x 268c?) were used to ensure that fingers 
remained uncorrelated across the box height. In sets 5 and 6 (the two lowest-diffusivity 
runs), lower finger heights allowed for shorter computational boxes (67c? x 67c? x 67c?) but 
a higher resolution (192 3 grid points) was used to resolve the finer-scale structures. Near 
Rq = 1, where the system dynamics are most turbulent, one additional simulation with twice 
the vertical resolution (192 or 384 points vertically) was run for each set. In all cases, the 
measured fluxes in the more fully-resolved run differed by no more than a few percent from 
the lower-resolution run, confirming that the spatial resolution used was sufficient to extract 
accurate fluxes. Sample snapshots of two of our simulations are shown in Figure [TJ As an 
impor tant note, we also fi nd that, except when R Q — > 1/r, the finger aspect ratio is close to 



unity (iDenissenkovl l2010h . 



We define non-dimensional heat and compositional fluxes through the Nusselt numbers 
Nut and Nu M , ratios of the total flux to the diffusive flux of the field considered. Equivalently, 
Nut, m — 1 represents the ratio of the turbulent diffusivity to the microscopic diffusivity. We 
measure Nut and Nu^ from each simulation once the system has settled into a statistically- 
steady state. The results are summarized in Figure [2j 

For ease of comparison between different simulation sets, we defined the rescaled density 
ratio as r = (Rq — l)/(r _1 — 1) so that the instability range is r G [0, 1] in all cases. We 
find that all simulation sets collapse onto a single universal profile for the turbulent fluxes 
Nut(?") — 1 and Nu^(r) — 1 as Pr and r decrease, profiles which can be expressed as 

Nu T (r)-l = r 3/2 VP^f(r), (12) 

Nu„(r)-1 = \ —g(r), (13) 
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where the universal functions f(r) and g(r) are adequately fitted to the data using 

g(r)J(r)~ae- br (l-ry. (14) 

For f(r) the data suggests that a = 264, b = 4.7, c = 1.1, and for g(r) we find a = 101, 
b = 3.6, c = 1.1. Figure [2] shows these fits together with the rescaled Nusselt numbers. The 
physical interpretation of these scalings with r and Pr is the following: for u, <C Kt, the 
instability is driven by the compositional field and the dependence of compositional transport 
on kt must drop out. The turbulent compositional diffusivity = (Nu /t — 1)/^ is then 
proportional to the geometric mean of v and k^. The scaling for the temperature Nusselt 
number then follows by as suming that the nondimensional turbulent flux ratio 7 =< wT > 
I < wp, > is always 0(1) (jRadkjbooj ). 



Although these simulations were conducted at parameter values a few orders of mag- 
nitude higher than those appropriate for stellar interiors, we find that they follow universal 
asymptotic laws. The simplicity of the system considered (e.g. no boundaries) also sug- 
gests that these laws can reliably be extrapolated to the astrophysical regim e (as long as 
Pr is of order r). From this we draw two critical observations. As found by iDenissenkov 
( 120101 ). the contribution of fingering convection to heat transport at stellar parameter values 



(Pr ~ 10~ 6 , r ~ 10 -6 ) is entirely negligible. It therefore does not affect the thermal structure 
of the object in any way. Compositional transport, on the other hand, can be significantly 
enhanced above molecular diffusion, by about two orders of magnitude depending on the 
ratio Pr/r = v/k^. 



W e now compare our flux laws with the parameterization of lUlrichl (Il972l ) and lKippenhahn et al 
(Il980l ). in which Nu^ - 1 = C P (L/W) 2 /tR* , where L/W is the aspect ratio of the fingers 



(leng th L divided by width W) and C p i s a model constant (C p = 1 for iKippenhahn et al. 



19801 . C p = 87r 2 /3 ~ 26 for lUlrichl Il972l) . Figure [ 3] com pares as a function of density 
ratio, at the fluid parameters used by iDenissenkovl (120101 ). Pr = 4 x 10~ 6 and r = 2 x 10~ 6 . 
Given that the numerically-deter mined finger a spect ratio is close to one for all simulations 
as r <^ 1, Pr <^ 1, we find that the lUlrichl (jl972j) prescription fares comparatively better than 
Kippenhahn et al.l (119801 ) . However, it still vastly overestimates mixing both near R = 1 
(where the system is closest to Ledoux instability) and Ro = 1/r, close to linear stability. 
Crucially, our results show that the diffusivity obtained is not sufficient t o explain RGB abun- 



danc e observations for w hich a value of C P (L/W) 2 ~ 1000 is needed (jCharbonnel fc Zahn 



20071 : IDenissenkovl l2010h 
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3.2. Layer formation 



The small-scale flux laws presented above are accurate measurements of transport by 
homogeneous fingering convection. However, they may not appropriately describe trans- 
port in the presence of the kind of large-scale thermo compositional stairc ases which are 



known to form in oceanic thermohaline convecti on (e.g . ISchmitt et al.ll2005l ). The origin of 
such staircases has recently been discovered by iRadko (£3), who demonstrated the exis- 
tence of a positive feedback mechanism that can, in some circumstances, drive the growth 
of horizontally-invariant perturbations in the temperature and salinity fields: large-scale 
variations of the density ratio lead to convergences and divergences of fingering fluxes that 
reinforce the original perturbation, producing a growing disturbance that ultimately over- 
turns into regular "steps." This mechanism relies fundamentally on the variation of the flux 
ratio 7 = F T /F S (where Ft is the turbulent heat flux and F$ is the turbulent salt flux) with 
density ratio, and was therefore called the 7-instability. The 7-instability mechanism accu- 
rately predicts the growth rate of large-scale perturbatio ns, and subsequent overturning into 



layers, in both two-dim ensional simulations dRadkoJ 2003h . and thr ee-dimensional simulations 
of staircase formation (ITraxler et al.ll2010l ; IStellmach et al.ll2010l ). The instability requires 
that 7 decreases as Rq increases, a condition which must be experimentally determined from 
the fluxes for each set of fluid parameters Pr and r. 

In Radko's original formulation for oceanic staircases, the systems of interest were dom- 
inated by turbulent fluxes. Diffusive contributions to the total fluxes were neglected for 
convenience and simplicity. However, as shown in Section 13. 1[ turbulent heat transport is 
negligible in the astrophysical case, so the 7-instability theory is re-derived here including 
all diffusive terms for accuracy and completeness. We begin with the equations of Section [2] 
and average them over many fingers: 



1 fdu „ 
— — + u • Vu 

Pr V dt 

dT „ m 
— + w + u • VT 

at 



1 



- Vp + (T — /j,)e z + V 2 u - —V • R 



Pr 



dfi 1 

m + R7 w + u 



V/i 



V 2 T - V • Fn 



rV> - V • F 



(15) 

(16) 
(17) 



where u, T, and \i now represent averaged large-scale fields. On the right-hand sides appear 
the usual Reynolds stress term, Rij = UjQj, and the turbulent fluxes = uT, F M = ujl. 



We proceed as in earlier formulations ( lRadkoll2003t ITraxler et al.ll2010l ) by making sev- 
eral simplifying assumptions, the first being that Reynolds stresses are small enough to 
neglect and the second being that the turbulent fluxes are dominated by their vertical com- 
ponent, so that Ft ~ -Fr e z ; F„ w Fne z . The 7-instability involves only the temperature 
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and compositional fields, so we consider only zero-velocity perturbations and discard the 
momentum equation. We then define the total fluxes as well as their ratio: 



77lt0t 

b T 

Tjitot 

7 



tot 



F T -{l + dT/dz), 
F„-T(l/Ro + dfx/dz), 
F^ ot 



(18) 
(19) 

(20) 



so that the thermal and compositional Nusselt numbers are: 

jptot 

Nu T = — ^ - ; , Nu 



77tot 



— (1 + dT/dz) 



-T(i?o + dfi/dz) 



The final assumption is that Nur, Nu M , and 7 depend only on the local value of the density 
ratio R p . Note that R p 7^ Rq and may vary with z. We now linearize equations f)16p and 
ifPTj) around a state of homogeneous turbulent convection in which T = + T", \i = + //, 
and R p = Rq + R' . For example, linearizing the density ratio R p , we have: 

aT 0z (l + dT'/dz) 



Rn 



The temperature equation yields 
dT' dF^ ot 



R (l + dT'/dz - Rodfi'/dz). 



(21) 



Of 



dz 
<9Nu T 



dz 
dNu T 



+ Nu- 



d 2 r 



0R P 



Ro 



dRp 

dz 



dz 2 

+ m T (Ro 



OR, 



Rq 



Ro 



d 2 r 

dz 2 



Ro 



) d 2 r 

dz 2 

ay 

dz 2 



+ Nu T ( J R 



d 2 r 

dz 2 



and similarly for equation ( TT7|) . Assuming normal modes of the form T',fi' ~ e 
tions ([221 and the equivalent linearization of (fT7|) combine into the quadratic 



ikz+Xt 



A 2 + Afc 2 



Nu (l - AM + A 2 1 



T tot 



- fc 4 AiNuoi2o = 0, 



(22) 



equa- 



(23) 



where we use the following notation for simplicity: 



A x = Rt 



d(l/-Y tot ) 



dR a 



An 



Nu = Nu T (i? ), 



Ro 



dNu T 



(1R 



Ro 



7o° t = 7 tot (^o)- 
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Note that this is identical to the quadratic of iRadkol (120031 ) replacing 7 by 7 



tot 



As discussed by IRadkol (120031 ). a sufficient condition for the instability is that 7 tot is 
a decreasing function of density ratio. Indeed, when this is the case, the coefficient A\ 
is positive, and the expression (123]) has two real roots, one positive and one negative. In 
order to determine whether the 7-instability occurs and therefore whether staircases may 
spontaneously form at low Pr, low r, we simply need to calculate 7 tot from our previous 
measurements. The results are shown in Figure H] We find that the flux ratio always 
increases with density ratio. Crucially, this suggests that thermocompositional layers are 
not expected to form spontaneously in astrophysical fingering convection. In the absence of 
such "staircases," the fluxes are accurately supplied by our flux laws [12] and | 



Conclusion 



Our results can be summarized in a few key points. Using high-performance three- 
dimensional numerical simulations we find that turbulent compositional transport by finger- 
ing convection follows the simple law: 



Dfj, = k m (Nu m - 1) = 101 v /^z7e" 



-3.67V 



1 1.1 



(24) 



where r = (Rq — lVf r -1 — 1). Meanw hile, the turbulent heat flux is negligibly small at 
stellar parameters (see lDenissenkovll20 1 Or ) . We also find that large-scale thermocompositional 
staircases are not expected to form spontaneously for Pr , r <C 1, so that the turbulent 
diffusivity proposed above may and should be used as given to model astrophysical fingering 
convection. 

Applying these findings to the problem of RGB stars, we concur with Denissenkov's 
conclusion that mixing by fingering convection alone cannot explain the observed RGB 
abundances, and that additional mechanism s should be investigated instead (e.g. gyroscopic 
pumping, see lGaraud fc Bodenheimerll2010l ). A second exa mple of applic ation to the plane- 



tary pollution problem is presented in a companion paper (lGaraudll2010l ) 
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Figured] was rendered using VAPOR (jClyne et al 
of the National Center for Atmospheric Research. 
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Clyne fc Rast 



20051 ). a product 



- 10 - 



REFERENCES 

Baines, P., & Gill, A. 1969, J. Fluid Mech., 37 

Calzavarini, E., Doering, C. R., Gibbon, J. D., Lohse, D., Tanabe, A., & Toschi, F. 2006, 
Phys. Rev. E, 73, 035301 

Charbonnel, C, & Zahn, J. 2007, Astron. Astrophys., 467 

Clyne, J., Mininni, P., Norton, A., & Rast, M. 2007, New J. Phys, 9 

Clyne, J., & Rast, M. 2005, in Proceedings of Visualization and Data Analysis 2005 

Denissenkov, P. A. 2010, ArXiv e-prints 

Garaud, P. 2010, submitted to ApJ Letters 

Garaud, P., & Bodenheimer, P. 2010, ApJ, 719, 313 

Kippenhahn, R., Ruschenplatt, G., & Thomas, H. 1980, A&A, 91, 175 

Kunze, E. 2003, Prog. Oceanogr., 56, 399 

Radko, T. 2003, J. Fluid Mech., 497, 365 

Schmitt, R. 1979, Deep-Sea Res., 26A, 23 

-. 1983, Phys. Fluids, 26, 2373 

Schmitt, R., Ledwell, J., Montgomery, E., Polzin, K., & Toole, J. 2005, Science, 308, 685 
Stancliffe, R. J. 2010, MNRAS, 403, 505 

Stellmach, S., Traxler, A., Garaud, P., Brummell, N., & Radko, T. 2010, ArXiv e-prints 
Stern, M. 1960, Tellus, 12, 172 

Stern, M., Radko, T., & Simeonov, J. 2001, J. Mar. Res., 59, 355 

Traxler, A., Stellmach, S., Garaud, P., Radko, T., & Brummell, N. 2010, ArXiv e-prints 
Ulrich, R. K. 1972, ApJ, 172, 165 
Vauclair, S. 2004, Astrophys. J., 605, 874 
Veronis, G. 2007, Deep-Sea Res. I, 54 



This preprint was prepared with the A AS IATgX macros v5.2. 



- 11 - 




Fig. 1. — Compositional perturbations, in units of the compositional contrast across the 
height of the domain, of two simulations at Pr = r = 1/10, at density ratio Rq = 1.45 (left) 
and Ro = 9.1 (right). In general, fingering convection is more laminar, and the fingers are 
taller as Rq — >■ 1/r (as the background stratification becomes more stable). However, as Pr 
and r decrease, secondary instabilities enforce shorter fingers even at high density ratio. 
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Fig. 2. — Turbulent heat and compositional fluxes as a function of rescaled density ratio r. 
Dividing the turbulent heat flux by r 3 / 2 Pr 1//2 and compositional flux by y / Pr/V reveals the 
universal scaling laws (fl2"|) and (|T3|) . Also shown are f(r) and g(r) fitted to the data (solid 
lines) . 
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Fig. 3. — Comparison between o ur empirically-dete r mined com positional tur bulent diffusiv 
ity and the parameterizations of iKippenhahn et al.l ( 1980 ) and lUlrichl ( 1972) (where the co- 



efficie nt C P (L/W) 2 serves as a free parameter). Since we find L/W ~ 1, IKippenhahn et al. 



(Il980l ) unde restima t es D ,,, for most of the Rq range and overestimates it as Rq — > 1 and 



Rq —t- 1/r. lUlrichl ( 119721 ) fares slightly better at intermediate R . Our results show that 



turbulent mixing by fingering convection is not sufficient to explain observations, for which 
C P (L/W) 2 ~ 1000 is required (see text for detail). 
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Fig. 4. — The total flux ratio 7 tot = i^ ot /F* ot for various systems. As the density stratifi- 
cation becomes more stable (r — > 1), the fluxes of heat and composition are dominated by 
diffusive contributions and their ratio approaches the limiting value Rq/t. The flux ratio 
increases monotonically with r, indicating that staircases will not form spontaneously in this 
parameter regime (Pr < l,r< 1). 
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Table 1: Summary of governing parameters for all simulation sets. 



Set 


Pr 


T 


R 


1 


1/3 


1/3 


1.025, 1.075, 1.1,1.1°, 1.1/, 








1.125,1.2,1.3,1.5,2,2.4, 








2.8,2.9 a - a 


2 


1/3 


1/10 


1.5, 1.5", 4, 7,9.5 a ' a 


3 


1/10 


1/3 


1.1,1.1°, 1.7, 2.3, 2.8 a ' Q , 








2.9 a ' a 


4 


1/10 


1/10 


1.1,1.45, 1.45°, 2, 2.8, 3, 








3.3,5,7,9.1 a - a 


5 


1/10 


1/30 


ll 6 ' 7 , ll 6 ' 5 , 20 6 ' 7 


6 


1/30 


1/10 


46,7 7 b >7 



Note. — Relevant parameters include values of Pr and r, of the density ratio i? € [1,1/r], domain 
size, and resolution. Unless otherwise specified (see footnotes), simulations were run at domain sizes of 
67 x 67 x 107.2d and resolution of 96 3 grid points. 

a 83.75(2 x 83.75(2 x 268(2 

fc 67(2 x 67(2 x 67(2 

"96 x 96 x 192 

^192 x 192 x 384 

n92 x 192 x 192 

5 384 x 384 x 384 



